Local Defect in Metallic Quantum Critical Systems 
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We present a theory of a single point, line or plane defect coupling to the square of the order 
parameter in a metallic system near a quantum critical point at or above its upper critical dimension. 
At criticality, a spin droplet is nucleated around the defect with droplet core size determined by the 
strength of the defect potential. Outside the core a universal slowly decaying tail of the droplet is 
found, leading to many dissipative channels coupling to the droplet and to a complete suppression 
of quantum tunneling. We propose an NMR experiment to measure the impurity-induced changes 
in the local spin susceptibility. 
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The behavior of 'droplets' of local order in a non- 
ordered background is an issue of wide relevance in con- 
densed matter physics. One particularly interesting sub- 
class of problems concerns 'droplets' induced by defects 
in nearly critical systems. A long-standing problem in 
heavy fermion physics concerns the very small magnetic 
moments which have been observed in several materi- 
als [0J^] and may be related to grain boundaries and 
other structural defects |^-||. In a colossal magnetore- 
sistance material, magnetic order was observed to be en- 
hanced near grain boundaries. Q A related issue is the 
magnetism induced in high temperature superconductors 
by apparently non-magnetic substituents such as Zn , 
which have been interpreted |^,^| as spin droplets induced 
in a nearly critical system (although other interpretations 
exist also [|l0)). Nucleation of regions of charge density 
wave order around defect sites on the surface of a 'corre- 
lated' material was reported by pi. 'Quantum Griffiths' 
effects and 'Kondo disorder' are presently of intense in- 
terest. ]l^l5|. The problem bears on the fundamental 
issue of the Kondo effect near a quantum critical point 
Jl6| . Finally, recording of information involves the polar- 
ization of small domains, whose long time dynamics and 
stability are of great importance. 

This Letter presents the theory of the local polariza- 
tion ('droplet') induced by a single defect in an otherwise 
nondisordered system which is near a quantum critical 
point at or above it upper critical dimension, d u . These 
restrictions allow a controlled theoretical treatment and 
apply to a wide range of systems including metallic mag- 
nets in dimensions d = 2, 3 [ |l7| , |l8| and "quantum para- 
electric" (i.e. nearly ferroelectric) systems in d = 3 |L9) . 
We study defects which couple to the square of the or- 
der parameter, i.e. change the 'local T c \ and thus create 
small regions where order is more favored than in the 
pure system. We address three questions: under which 
circumstances does the defect create a 'droplet', a re- 
gion about the defect in which the order parameter is 
non- vanishing (at least on short time scales)? What is 
the size and general properties of the droplet? What are 
the relevant fluctuations? Our work is complementary to 



that of Vojta and Sachdev p3| , who studied a linear cou- 
pling of the defect to the order parameter in a quantum 
critical system below d u . It is also related to the work 
of Castro-Neto and Jones [jl3 14 1, but differs in several 
important aspects. 

Our starting point is a quantum Ginzburg-Landau ac- 
tion for an order parameter field <f>. After obtaining its 
mean field solution we consider fluctuation corrections, 
which are tractable because the dimension of the system 
is above d u . Halperin and Varma used a similar approach 
for a classical system pl| . Our action, in conveniently 
scaled variables, is, in space and imaginary time: 

S = S 4y*+lf " dT J d d x {{V {*) + k 2 ) 0(x,r) 2 



+ (V0(x,r)) 2 + ^(x,r) 4 
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Here, k determines the distance of the bulk system to the 
critical point. We measure lengths in units of the bare 
correlation length of the problem (typically of the order 
of a lattice constant) and energies in terms of the con- 
densation energy Eq obtained by evaluating the static, 
spatially uniform free energy with n = 1. In the fol- 
lowing, we consider only symmetrical defects, which are 
characterized by a dimensionality dd (i.e., the number 
of dimensions along which the defect potential, V, re- 
mains constant), a length scale a (expected to be of 
the order of a lattice constant) over which V decays in 
the D = d — dd transverse directions, and a dimension- 
less strength v = — j d D rV(r) (v > corresponds to 
a local tendency towards order). Here, r refers to the 
D = d — dd transverse components of the d-dimensional 
vector x, where dd = 0, 1,2 corresponds to a point, line 
and plane defect, respectively. 

The dynamic term <Sd yn takes the general form 
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where the coefficients c and T depend on whether the 
system is overdamped (metallic case) or not, on the sym- 
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metry of the order parameter, and on whether it is con- 
served. Examples include: (1) the undamped Ising an- 
tiferromagnet, with = and c = const; (2) the un- 
damped Ising ferromagnet with a conserved order param- 
eter, r^ 1 = and c ~ l/q; (3) the metallic (overdampcd) 
antiferromagnet T = const; (4) the metallic Ising 

ferromagnet r q ~ q. The form of <Sd yn combined with 
the static part of the free energy defines a (mean-field) 
dynamic exponent z which is z = 1,2,2,3 respectively 
for the cases listed above. The effective dimensionality 
of the quantum phase transition problem defined by S is 
d e g = d + z and we restrict to d e g > 4. 

We now sketch the essential features of the mean field 
solution (details will be given elsewhere p2|). We focus 
only on the transverse dimensions and assume V(r > 
a) = so that for r > a the mean field equation is 



V 2 o + n 2 (t> + 4>l = 0. 



(3) 



For < k <C 1, the solution is of the form <fio(r) = 
r^ 1 f(r /tq, ktq) where / is dimensionless. The length 
scale ro is determined by connecting the solution of 
Eq.(||) (i.e., for r > a) to the one for r < a and thus 
depends on the defect strength, v. For nr > 1 the </> 3 
term may be neglected at all r and the solution is the 
familiar exponentially decaying one. In the more inter- 
esting case Kro < 1 and D < 3, the behavior at r < k" 1 
is controlled by the nonlinearity and the scale ro defines 
the size of the droplet. We find (the go are constants): 
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and in D = 3 up to logarithms, the 
ro < r < k~ l behavior of c/)q (r) is independent of the 
short length scale physics encoded in ro- In all dimen- 
sions J d D r(po(r) diverges at criticality and in D > 2 
J d D r(j>l(r) diverges. Thus many physical properties are 
dominated by the exponential tail of <fio{r). Fig. [j]a 
shows a schematic picture of the droplet amplitude. 

The scale ro may be estimated by substituting our re- 
sults for cf>Q (r) into Eq.(Q) and minimizing with respect 
to r . This also determines the binding energy, -Ebind, 
of the droplet which forms at temperatures T < -EWd. 
Alternative approaches including scaling, exact solutions 
(D = 1), and numerics give identical results |p2| . In all 
cases the energetics are dominated by the droplet 'core' 
region r < rg. At criticality, an arbitrarily weak po- 
tential induces a droplet in D = 1,2 but in D = 3 
a critical strength is required. In D = 1 ro 
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-Eqv a while in D = 2 ro 



and 



-Ebind ~ —E ve 2 I'" In D = 3 a critical value v* 
is required for droplet formation, and for v > v" 
is of order a while £ b ind ~ -E ^r(^ - l)a~ 2 . 
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FIG. 1. (a) Schematic r-dependence of the droplet am- 
plitude, 4>o(r). In D — 2, <f>o(r) = const, for r < ro up 
to logarithmic corrections, (b) Dependence of even parity 
wave-functions for k — 0.5 and k = 1.5 on distance from the 
defect. Inset: near-defect (r w ro) region for k = 0.5. 

droplet size cannot exceed which yields an estimate 
for the critical potential v*(k) v > v* = 2k (D = 1) and 
v > v* = — 27r/log(Ka) (D = 2). The droplet magneti- 
zation M c i in the ferromagnetic case is given by the in- 
tegral of 4>q (r) which diverges as n l ~ D (logarithmically 
in D = 1) as criticality is approached. In an antifer- 
romagnetic system with characteristic wave vector Q, 
M d ~ |Q| which for D > 1 is a number of order 
unity even at criticality. We show below that these mean 
field results are not significantly affected by fluctuations, 
at least in metallic systems. 

Gaussian fluctuations may be treated by expanding to 
quadratic order about the mean field solution, leading to 
the action: 
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where L = k 2 + V (r) + 3</>q ( r ) — V 2 .. L has only posi- 
tive energy eigenfunctions. At criticality all of these are 
extended but if k then for v in a small range above 
v*, a bound state of energy < E < k 2 may occur. The 
form of the potential (weak slowly varying repulsion with 
attractive center) leads to non-monotonic wave-functions 
with an upwards cusp at the defect scale a, a decrease 
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with distance in the range between a and Tq, and then 
(for extended states) an increase back to the unit ampli- 
tude of a propagating plane wave. This r— dependence is 
shown for two even parity wave-functions of the Gaussian 
fluctuations in Fig. 0b. The effect of Gaussian fluctua- 
tions on the droplet size and shape may be computed 
in terms of the difference between the eigenfunctions of 
L and of k 2 — V 2 and are found (2^] not to change the 
long distance behavior of the droplet (as expected from 
conventional RG arguments |l7| , [l8|] ) . 

A more important class of fluctuations changes the ori- 
entation of the droplet. In a compact {6,4 — 0) droplet 
these are rotations and 'instanton' processes in which the 
droplet collapses and re-forms. In extended (dd > 0) 
droplets the important processes are motion of domain 
walls. Here, we sketch the results for a compact droplet; 
details and an analysis of the moving domain wall case 
will be given elsewhere ^2|. The dynamics of an iso- 
lated (not embedded in a critical system) droplet have 
been previously studied the new feature here is the 
overdamped dynamics (in the metallic case). 

We first consider a droplet with Ising symmetry, in 
which case the instanton, i.e., the collapse of the droplet, 
is the important fluctuation process. To estimate the 
action we substitute the ansatz 4>(r,r) = 4>o(r)r](T) into 
Eq.(|l]) and retain leading time derivatives. We find that 
the action corresponding to 2N instantons is: 
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S 2 n = ^ £ log( W -ife)(-l) i+ ' + 2JV ^ + 

+ J du J dvlog(l + y*(u-v) 2 ), (5) 



with C - 15£o/(4£bin d ), m = £ 2 E q Cq 2 <Mq) 2 and 
7 = E X^qTqVc^q) 2 - The terms proportional to 2N 
are the action of a single instanton and yo (determined 
by minimizing the second and third term of S^n) is the 
duration of an instanton in units of E^ ind . Note that both 
m and 7 diverge as k — > and that 7 = in undamped 
models. In the weak dissipation (7 <C 1) limit the stan- 
dard macroscopic quantum tunneling analysis [Q leads 
to S(yo) — 8^m/( + 6m'y( + (j 2 ) so the 'bare' droplet 
tunneling rate is ~ Eoe~ s< ' v °' > and vanishes as criticality 
is approached. Instanton-instanton interaction effects, 
which arise from the first term in S2N are handled via 
a perturbative renormalization group treatment and if 7 
is less than a critical value ~ 1 dissipative effects reduce 
the tunneling rate but not to zero. 

In a metallic system near criticality, 7 ^> 1, and the 
conventional analysis does not apply. Our detailed re- 
sults depend on the ratio m/7. If 771/7 ^ 1 (f er ~ 
romagnetic case) then minimization leads to S(yo) = 
27 [log(m/27) + 1] while for 771/7 « 1 w find S(yo) = 
1/3 

3m(7/6m) ' . In cither case, dissipation strongly sup- 
presses the bare tunneling rate, and the large value of 
7 puts the action on the localized side of the Caldeira- 



Leggett phase boundary, implying that tunneling pro- 
cesses are completely suppressed on long time scales. 

Our treatment closely parallels Hamann's formulation 
pB| of the Kondo dynamics of a single spin in a non- 
critical metal. Hamann found S(yo) = ln(l/JiVo) (J is 
the Kondo coupling and Nq is the Fermi surface density 
of states) and 7 = (1 — JNq) 2 . The crucial difference 
is that in our problem the large size of the droplet al- 
lows many dissipative channels to couple to it, leading 
to much stronger dissipative effects. Castro-Neto and 
Jones jl3| argued that the tunneling of a droplet could 
be mapped onto a single-channel Kondo problem; in our 
model this is not the case. A subsequent paper jl4j con- 
sidered a droplet consisting of a large number of elemen- 
tary 5 = 1/2 spins locked together by some magnetic 
interaction in a nearly critical system neglecting the \/r 
tail of the droplet and studying the effect of dissipation on 
the bare tunneling rate. Their specific results therefore 
differ from ours. Their crucial findings were that for anti- 
ferromagnetic systems droplets as large as I0 3 spins could 
tunnel and that dissipative effects only become impor- 
tant below an exponentially small scale, leaving a wide 
regime where quantum Griffiths behavior occurs, whereas 
we find that near criticality dissipation always dominates 
and quantum effects are suppressed. 

In droplets with XY or higher symmetry, rotational 
fluctuations must be considered. We integrate out the 
conduction electrons of an action of the type discussed 
in p6[ and assume that <f> is characterized by an ampli- 
tude <fio (r), obtained by solving Eq.(^), and a direction 
n(r) = (cos(0(t)), sm(#(r)), 0) specified by an angle 9. 
Expanding in the angular variables and retaining leading 
time derivatives gives 
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J dTidT 2 d T n(T 1 ) ■ <9 t ii(t 2 ) In ' 11 1 
J dr (id r n x n-z + h z {r)) 2 
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with M 2 = J d D rcj) 2 (v) h(r) = 2f d d r<j) {r)54> z (r,T)lM 2 
and ^ xy = (27r)~ d J d d q \(f>o(q) 2 \ /T q . We thus obtain the 
action expected for a rotor with a large moment of in- 
tertia (M) in a dissipative environment, precessing in an 
effective magnetic field caused by the background spin 
fluctuations. Ref. [g7| indicates that in the large dissipa- 
tion limit the subtleties associated with spin quantization 
may be neglected. A straightforward variational estimate 
then yields sub-diffusive long-time behavior, correspond- 
ing to a divergent susceptibility. 

The presence and fluctuations of the droplet are in 
principle observable via NMR measurements of the spin 
lattice relaxation rate, T± (r) ~ TImx(r, uS)/uj and local 
Knight shift. There are two different contributions: from 
changes, due to the droplet, in the extended 'Gaussian' 
spin fluctuations, and from the presence and tunneling of 
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FIG. 2. Dependence of the Gaussian fluctuation contribu- 
tion to the NMR relaxation rate on distance from the defect, 
calculated for z=2, d=2 and D=l and two distances from 
criticality. Inset: expanded view of the near-defect (r « ra) 
region. 



the droplet itself. The two contributions have very dif- 
ferent timescales and position dependences. The Gaus- 
sian fluctuations have characteristic frequency lo ~ n z 
and give rise to the usual divergence of the bulk relax- 
ation rate as criticality is approached. This contribu- 
tion to 1/7]T is suppressed near the defect (as shown 
for example in Fig. 2) essentially because the droplet 
reduces the amplitude of the low energy wave functions 
in the near-defect region (cf Fig. lb). The droplet tun- 
nelling processes provide a contribution to x <*0 which 
is proportional to the square of the droplet amplitude but 
varies on a much slower timescale, which moreover van- 
ishes exponentially as criticality is approached, so that 
these processes drop out of the NMR frequency window, 
appearing instead as a broadening oc 1/ro of the NMR 
spectrum. Further details will be given elsewhere [p2f . 

In summary, we have presented a theory of a single 
defect in a quantum critical system at or above its upper 
critical dimension. A crucial property is the 1/r 'tail' of 
the droplet extending into the surrounding medium. The 
ease with which line and plane defects induce regions of 
local order may be relevant to the small moments ob- 
served in heavy fermion systems [fj]j|]. Our finding that 
near criticality in a metallic system droplets behave in 
an essentially classical manner leaves no significant pa- 
rameter regime in which the quantum Griffiths behavior 
discussed in |jl3 14 exists. The strong dimensionality de- 
pendence of our results has implications for the general 
issue of Griffiths behavior near quantum criticality. Fi- 
nally, our results bear on the fundamental question of the 
Kondo effect near a quantum critical point A sin- 

gle spin in a nearly critical system will similarly induce 
a large droplet, which we believe will be prevented from 
tunneling by dissipative effects . 
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